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1 Introduction 



In the Hamiltonian formalism one can calculate the spectrum directly jl], || . In the 
Coulomb gauge, the color form of Gauss's law, which would appear to be essential for 
confinement, is satisfied exactly. Moreover an obvious confinement mechanism suggests 
itself in the Coulomb gauge namely, a long-range instantaneous color-Coulomb potential. 

This raises the question, can the Coulomb gauge be regularized? Certainly, on the 
lattice any configuration may be fixed to the Coulomb gauge. The lattice Coulomb hamil- 
tonian which will be described in this lecture is an implementation of the Coulomb-gauge 
dynamics with lattice regularization. As to whether the Coulomb gauge is perturbatively 
renormalizable in the continuum theory, the question is under active investigation |5|. || . 
The answer appears to be yes. 

2 Classical Yang-Mills Equations 

Classical Yang-Mills theory is designed to be invariant under the group of local gauge 
transformations g(x). We shall consider the local SU(N) theory, g(x) G SU(N). The 
quark field transforms covariantly, if)(x) — > 9 ip(x) = g(x)ip(x). The gauge-covariant 
derivative D = D(A), defined by {D lJi ip){x) = [<9 M + t a A < t{x)\'ip{x) also transforms covari- 
antly, (D ft i/))(x) — > g(x)(D fl i/))(x), provided that the connection A transforms according 
to 

A^x) ^A^x) = g{x)d lx g\x) + g{x)A tl {x)g\x). (1) 

Here the anti-hermitian matrices t a form a representation of the Lie algebra of SU(N), 
[£ a ,t 6 ] = f abc t c , and A^ = t a A a ^. The local gauge principle states that A and 9 A are 
physically identical A ~ 9 A The physical configuration space Q = A/Q is the space of 
connections A = {A} modulo the local gauge transformations Q = {g}. This space is 
non-trivial because the local gauge transformation is non-linear. 
The Yang-Mills field tensor, 

= d^A v - d u A^ + [Ap, A v \ = eF^, (2) 

transforms homogeneously, 

F^x) ^F^x) = g{x)F llv {x)g\x). (3) 

The Yang-Mills action 

Sym = 9 o 2 J d'x(l/A){-g^g^)F; u F^ (4) 
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is gauge-invariant and yields the field equations D^F^ = d^F^ v + [A^, F^ u ] = 0. Here 
the metric tensor is defined by g^ v = diag(— 1, 1, 1, 1), and g is the coupling constant. 
(If quarks were present we would have D^F^ U = — j qu ,u-) 

The color-electric and magnetic fields are defined by Ei = F i0 = —Ai + DiA and 
B, t = F jk for k cyclic. 

3 Gauss's law and Color- Coulomb Potenial 

The time derivative of A nowhere appears in the preceding equations, and the field equa- 
tion which corresponds to variation with respect to A is the time-independent constraint 
which constitutes the color Gauss's law, 

DiEi = ViEi + Ai x Ei = (or = p qu ), (5) 

where (A x E) a = f abc A b E c . This law may be used to fix the longitudinal part of the 
color-electric field. 
We write 

Ei = Ef - v t n. (6) 

where ViEf = 0, and Q is the color- Coulomb potential. We have 

DiEi = -DiViQ + Ai x Ef = p qu . (7) 

The three-dimensional Faddeev-Popov operator is defined by 

M = M{A) = -Di{A)Vi (8) 

Gauss's law fixes the color-Coulomb potential 

n(x) = (M- l p)(x) = J d 3 yM-\x,y,A)p(y). (9) 

This equation, which appears to be neither Lorentz nor gauge non-invariant, in fact holds 
in every Lorentz frame and every gauge. Here p is a color-charge density 

p = -AxE^ + p qu (10) 

which does not contain the complete color-charge density of the gluons, p g \ = —A x E, 
but only the part that comes from E tT . 

The color-Coulomb potential is propagated by the Green's function 

M- 1 = i-DiVi)' 1 = (-V 2 -Ax V,)- 1 (11) 



2 



of the three-dimensional Faddeev- Popov operator. The famous anti-screening property of 
non-Abelian gauge theory arises from this circumstance. For whereas (—V 2 ) is a positive 
operator, the second term (—Ai x V,) may be positive or negative. In the latter case 
there is a cancellation in the denominator which produces an enhanced color-Coulomb 
potential. 



4 Classical canonical formalism 

The Yang-Mills lagrangian is of the form 

L YM = J d 3 xC YM (12) 

where the Lagrangian density is 

£ YM = ^ 2 (l/2)(£ 2 -5 2 ). (13) 

The time derivatives are contained in = —Ai + D^Aq, and the canonical momentum 
densities are given by 

Tii = OC/dA l = -g^Ei 

7T = dC/dA = 0. (14) 

The classical Hamiltonian 

H = J (PxfcAi-C) 

H = jd 3 x[ (l/2)(s 2 vr 2 + g^B 2 ) - A D m ] (15) 

produces the equations of motion 

Ai = 5H/5TTi 7Tj = -g' 2 Et 
iii = -5H/5A, D Ei + e ljk DjB k = 0, (16) 

and the Gauss's law constraint 

= 5H/5A = (17) 
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5 Quantization in the Weyl gauge 



Before deriving the lattice Coulomb Hamiltonian, we shall derive the Christ-Lee Coulomb 
Hamiltonian for continuum gauge theory which it resembles. Christ and Lee start by 
canonical quantization in the Weyl gauge, A = 0. This is not a complete gauge-fixing 
and allows time-independent local gauge transforamtions g = g(x). In this gauge we 
have Ei = —A^ and the previous Hamiltonian reduces to 

H W = J d s x (l/2)( g 2 vf 2 + g 2 B 2 ). (18) 

This Hamiltonian generates the previous equations of motion in the Weyl gauge, but not 
the Gauss's law constraint D^i = 0. In the classical theory that is imposed as an initial 
condition that holds at t — 0, and which is preserved by the equations of motion. 

In the Christ-Lee method, one quantizes in the Weyl gauge by treating A and ir as 
Cartesian variables, with canoncal equal-time commutation relations 

[At(x),A b 0)] = K(x),n b M]=O (19) 



[<(£), A]{y)\ = -i5 ab Sij 5(x - y). (20) 

There is then no ordering problem in the Weyl Hamiltonian, nor in the definition of the 
left-hand side of Gauss's law 

G a (x) = (A^Or) = V*< + f abc A*>(xX(x). (21) 

Moreover with the above canonical commuation relations one observes that G a (x) is 
the generator of local time-independent gauge transformations, 

[G a (x)y i (y)]=i5(x-y)f* e *Z(y) (22) 

[G a (x), A\(y)\ = -iDf5(x - y), (23) 

under which 7r transforms homogeneously. The G a (x) satisfy the Lie algebra of the 
time-indpendent local gauge group, 

[G a (x), G\y)} = i5(x - y)f abc G c (y) (24) 

and moreover these local time-independent gauge transformations are a symmetry of the 
Hamiltonian [G a (x) , H\y] = 0. Therefore it is consistent to impose Gauss's law as a 
subsidiary condition which must be satisfied for a wave function ^ = ^f(A) to represent 
a physical state, 

G a (x)^ = 0. (25) 
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This condition is also the statement that the wave function is invariant under time- 
independent local gauge transformations. 

The inner product is defined by the functional integral 

($,*) = JdA $*(A)V{A), (26) 

and the Hamiltonian by Hw — T + V, where T and V correspond respectively to kinetic 
(electrical) and potential (magnetic) energies. Here 

V = fcfx (l/2)s - 2 [5?(*)] 2 (27) 

acts by multiplication, and we define T by its expression as a quadratic form 

($,T*) -JdA J d 3 x {gl/2) 5^{A)/5A a l {x) 5^>{A)/5A a i {x). (28) 

6 Minimal continuum Coulomb gauge 

In the Christ-Lee method 0], one solves the Gauss's law subsidiary condition explicitly 
to obtain the Hamiltonian in the physical subspace. In the original Christ-Lee paper, 
this is done by choosing coordinates in A-space adapted to the gauge orbit. A gauge 
orbit is a set of gauge-equivalent configurations A ~ 9 A. (Here we refer only to time- 
independent local gauge transformations g.) This may be done by writing A = 9 B, where 
B is transverse, V,-Bj = 0, and g runs over the set of all possible gauge transformations. 
Here the transverse configuration B labels the gauge orbit and g labels the point on 
the orbit. One way to address the problem of Gribov copies ||, is to restrict B to a 
subset A of transverse configurations V known as the fundamental modular region. For 
an alternative approach, see ||. 

A convenient may way to choose A is by means of a minimizing function 

F A {g) = || 9 A || 2 = J d 3 x\gV^ + gAg^\\ (29) 

The fundamental modular region A is defined to be the set of A that are an absolute 
minima of this function, 

A ee {A : || A \\ 2 < || 9 A || 2 forall g}. (30) 

With suitable assumptions on the space of configurations, one may show that such a 
minimum exists |H|. This definition (almost) uniquely fixes a single point on each gauge 
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orbit, to within global gauge transformations, because the absolute minimum of the 
functional is unique apart from "accidental" degeneracies. These are degenerate (equal) 
absolute minima of the minimizing functional. They constitute the boundary of A. The 
fundamental modular region is obtained from A by topologically identifying the boundary 
points that are gauge-equivalent degenerate absolute minima. (For a more extensive 
discussion of A see |Ll| .) 

We may easily derive some properties of the fundamental modular region. At a 
relative or absolute minimum the functional is stationary to first order and its second- 
order variation is positive. With g = exp(u;) = 1 + w + (l/2)u; 2 + ... an elementary 
calculation gives 

|| 9 A || 2 =|| A || 2 -2{WiA u u) + (1/2)(oj,M(A)oj) + ... , (31) 

where M(A) = —D^AjVi is the three-dimensional Faddeev-Popov operator. Since this 
must be positive for all a>, we obtain the Coulomb gauge condition VjAj = 0. In addition 
we find that M(A) is a positive operator, (u, M(A)u) > forall uj and for all A 6 A. 
These two conditions define the Gribov region Q which is the set of all relative and 
absolute minima, of which the fundamental modular region A is a proper subset, Acd. 
Note that because Faddeev-Popov operator M(A) is a positive operator for AG A, M(A) 
is invertible in the minimal Coulomb gauge and has non-negative eigenvalues. 

7 Continuum Coulomb- gauge Hamiltonian 

In the present lecture we shall not explicitly change variables from A to the adapted 
coordinates B and g, as was originally done by Christ and Lee j7|, but instead use the 
more efficient Faddeev-Popov "trick" to obtain i? C oui- 
The Faddeev-Popov identity reads 

1 = f dgd A (Vi 9 A)det[M( 9 A)], (32) 

where M(A) = — ViDi(A) is the Faddeev-Popov operator. Here 5x(ViAi) is the restric- 
tion of 5(ViAi) to the fundamental modular region A, namely 5a(ViA) = 5{V 'iAi)xA(A), 
where Xa(^4) — 1 for (transverse) A e A, and Xa(^) = otherwise. 

We require wave functions to be gauge invariant, \E f ( 9 A) = ^{A), and we wish to 
express the inner-product of such wave-functions ($, \I / ) as an integral over the parameters 
A tr G A that label the orbit space. For this purpose we insert the Faddeev-Popov identity 
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into the inner product and obtain 



/ dA &*{A)V{A) 




(33) 



Here N = J dg is the (infinite) volume of the gauge orbit which is an irrelevant normal- 
ization constant, and J A dA tr is the functional intregral over all transverse configurations 
in A. We have used the invariance of the measure dA under local gauge transformation, 
dA = dA' for A' = 9 Ai. These manipulations are somewhat formal in the continuum 
theory, but are well defined in lattice gauge theory. 

We next express the Hamiltonian Hyj = T + V on the reduced space of wave functions 
that depend on A tr . Because V(A) is gauge invariant, it acts on the reduced space as 
the operator of multiplication by V(A tr ). 

The expression for T as a quadratic form is particularly convenient for our purpose 
because it only involves first derivatives, and this allows us to apply Gauss's law. We 
have seen ( p2|) that vrf(x) = — i5/5A^(x) transforms homogeneously under gauge trans- 
formations, so for gauge-invariant wave functions the quantity £<3>* /5A^(x) 5*5/ /5A^(x) is 
also gauge- invariant. Consequently when we insert the Faddeev-Popov identity into the 
quadratic form that defines T, we obtain, as for the inner product, 



To proceed, we shall use Gauss's law Di5*5//5Ai = to express 5*5// 5 A^ as a derivative 
with respect to A tr . For this purpose we write Ei = Ef — Vifi, where Ei = i5/5Ai and 
Vj-Ej tr = 0. Here f2 is the color-Coulomb potential operator, and E a,tr = i5/5A^' tr . With 
DiEi = Vj-Ej + Ai x Ei, Gauss's law, reads 





x (g 2 /2) 5<5>*{A)/5A1{x) 5*5/ {A)/ 5 A* {x). 



(34) 



- AVifttf = p*5/. 



(35) 
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where we have introduced the color-charge density operator corresponding to the dynam- 
ical degrees of freedom 

p = -At x Ef, (36) 

As in the classical theory, we may solve for Qty by inverting the Faddeev-Popov operator 
M(A) = -Di(A)Vi. We have = M~ l p^, which gives 

Eft! = £f * - WiM-\A) (-Ai x Ef (37) 

Only derivatives with respect to A tr now appear in the expression for the quadratic form 
T, 

($,T^) = iV / dA tr det[M(A tr )] 

x J d 3 x {gl/2) (Ef(x)$)*(A tr ) (E?(x)^)(A tr ), (38) 

where E\I/ and E$ are defined by the preceding equation. The transverse and longitudinal 
parts of the electric field contribute separately to the kinetic energy, 

($,T$) = iV cL4 tr det[M(A tr )] J d 3 x {gl/2) 

x [ |£f' a (x)$| 2 + KViM-V)"^)^! 2 ]• (39) 

We have expressed T as a quadratic form on functions of transverse configurations, 
A tv G A that parametrize the orbit space, and V as the operator of multiplication by 
V(A tr ). These expressions define the Hamiltonian in the minimal Coulomb gauge, H con i = 
T + V, defined by the gauge choice Ai = Af 6 A. 



8 Kogut-Susskind lattice Hamiltonian 

In a Euclidean lattice theory the Hamilton H may be derived from the partition func- 
tion Z. The partition function is expressed in terms of the transfer matrix T by Z = 
liniTv^oo tvT N , where is the number of Euclidean "time" slices. One chooses an asym- 
metric lattice with lattice unit in the Euclidean time direction and lattice unit a in 
spatial directions. Then T = lim ao ^ exp(— a H). 

In |~T| the lattice Coulomb-gauge Hamiltonian was derived directly from this formula, 
where the configurations were fixed to the lattice Coulomb gauge. In the present lecture 
we present a much briefer and elementary derivation. We shall start with the known 
Kogut-Susskind Hamiltonian flT2| |, H^s, which is the lattice analog of the Weyl Hamilto- 



nian. We then fix the gauge, and solve the lattice form of the Gauss's law constraint to 
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obtain the lattice Coulomb-gauge Hamiltonian on the reduced space, just as we did in 
the continuum theory. 

We designate points on the three-dimensional periodic cubic lattice by three-vectors 
with integer components, which we designate x,y,... . Gauge transformations are site 
variables, g x G SU(N). The basic variables of the theory are the link variables U( xy ) G 
SU(N) defined for all links (xy) of the lattice, with U( yx ) = U^ xy y Local gauge transfor- 
mations are defined by U( xy ) -> 9 U( xy ) = g x U{xy)gl- 

The Uq — 1 gauge is the lattice analog of the Weyl Aq = gauge. In this gauge, the 
Hamiltonian obtained from Wilson's lattice action and the formula T = lim ao ^ exp(— aoH) 
is 

H KS = T + V (40) 

where 

t = J2 {glMJ 2 {xy) 

V = (2N/ag 2 ) ]T(1 - N^Re trU p ). (41) 
v 

Here a is the (spatial) lattice unit. Henceforth we set a — 1. In the last expression the 
sum extends over all plaquettes p of the (spatial) lattice, and U p is the product of the 
link variables U( xy ) on the links around the plaquette p. The potential energy operator 
is gauge invariant V( 9 U) = V(U), as are physical wave functions, ^( 9 U) = ^(U). 

The variables J( xy ) are electric field operators associated to each link of the lattice, 
defined as follows. Let the link variable U( xy ) be parametrized by a set of variables 0? xy \, 
then, for each link (xy), 

J {xy) ,a^iJ:(e {xy) )d/d9f xy) (42) 
is the Lie generator of the SU(N) group, that satisfies the commutation relations 

\J(xy),ai J(uv),b\ i$(xy),(uv) f J(xy),c- (^ > ) 

For the SU(2) group, T would be the Hamiltonian for a set of non-interacting tops, located 
on every link of the lattice. The potential energy operator V provides an interaction 
between tops on the same plaquette. 

Both T and V are gauge invariant, as is the inner product 

(¥,<&) = / HdU {xy) **(U)<S>(U). (44) 

(xy) 
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Here dU( xy ) is the Haar measure on each link, 



dU {xy) = det ^(xj/) II d() i 



(xy)i 



i a 



(45) 



a 



where il>(0( xy )) is the inverse of the matrix J(9( xy )), 



(46) 



9 Minimal lattice Coulomb gauge 



We wish to express the hamiltonian as on operator on the reduced space of gauge orbits. 
For this purpose we require a set of parameters that label the gauge orbits. A convenient 
way to identify the gauge orbits is to introduce a minimizing function, as in the continuum 
theory, 



Here we have introduced the notation U Xyi = U( xy ) for y = x + where e$ is a unit vector 
in the positive i-direction. The gauge-transformed link variable is given by ( 9 U) Xy i = 
9xUx,i9 x+ z, where we have written % = C{. The fundamental modular region A is defined 
to be the set of configurations U that are absolute minima of this function, 



With this gauge choice, the link variables U Xji are made as close to unity as possible, 
equitably, over the whole lattice. 

The advantage of this method is that we are assured that a minimizing configuration 
exists on each gauge orbit, because the minimizing function is defined on a compact 
space (a finite product of SU(N) group manifolds). An alternative procedure that is 
sometimes followed is to posit a gauge condition. For example with the exponential 
mapping U X j = exp^^j), one may posit the condition ^(Am — Ar-i) = ^ which is a 
lattice analog of the continuum Coulomb-gauge transversality condition. An obstacle to 
this procedure is that one does not know whether every gauge orbit intersects this gauge- 
fixing surface. On the other hand for the purpose of formal expansions the alternative 
procedure may be used. 

As in the continuum case, we observe that, at a relative or absolute minimum, the 
minimizing function is stationary to first order, and its second order variation is positive. 
In order to conveniently exploit these properties, we introduce a one-parameter subgroup 



F u (g) = J2{i-N- 1 RetT[^U) X;i ] }. 



(47) 



x,i 




(48) 
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of the local gauge group g x (s) = exp(su x ), where u x = t a uj% is anti-hermitian, uj x = —uj x . 
We write Fu(s) = F u (g(s)) and we have 

dFu{s)/ds = -N' 1 ^ tr[u x ( 9 U) X , - ( 9 U) x ,uj x+ - ] 

x,i 

= N- 1 ^ ^ tr{K + j - u x ](«U) Xii } (49) 

x,i 

= -(2JV)" 1 S trK[(^)„ - (W)t, -(V)^, + (*[/)t_, .] }. (50) 

At a relative or absolute minimum U we have dF[/(s)/ds| s= o = 0. To express this 
condition compactly, we introduce the link variables 

t C l A xi = A x<i = (1/2)(C/ Xj j — C/ Xj j) traceless- (51) 

We shall use these variables as coordinates to parametrize the SU(N) group, and we write 
U Xi i = U(A Xyi ). The manifold of the SU(N) group requires more than one coordinate 
patch, but we shall not trouble to account for this explicitly. The formulas of the last 
section hold, with 9 x i = A" v These coordinates agree with the exponential mapping, 
U Xj i = exp A Xj i, to second order, for with the above definition one may show 

U Xji = 1 + A X:i + (1/2)^) + (WyhrA 3 + 0(A 4 ), (52) 

so these variables approach the continuum connection A i {x) in the continuum limit. 
In terms of these variables, the condition dFu(s)/ds\ s= o = reads 

£lr[w*E(-^,i + 4.-^)1=0- (53) 

x i 

Since this holds for arbitrary u, we conclude that, at a minimum of the minimizing 
function, A satisfies 

E (Am - = 0- (54) 

i 

This is the lattice analog of the transversality condition that characterizes the continuum 
Coulomb gauge, and we call our gauge choice the "minimal lattice Coulomb gauge" . The 
lattice Coulomb-gauge condition is linear in the variables A Xti , and it may be solved by 
lattice fourier transform. 

To see the geometrical meaning of the gauge condition, and for future use, it is helpful 
to introduce some definitions. We define the lattice "gradient" by 

(Vw)«,i = (uj x+ --uj x ), (55) 
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It is a matrix that linearly maps site variables into link variables. The identity, 

^2tT[(u x+i - uj x )A x4 ] = ^2tT[u x ^2(-A Xii + A x _ tti )}, (56) 

x,i x i 

holds for all lj and A. We write it 

(Vw, A) = (u, V*A) (57) 

which defines the dual V* of the lattice gradient. It is a matrix that linearly maps link 
variables into site variables, that is given explicitly by 

(V*A) X = -J2(A X ,*-A x _ ll) ), (58) 

i 

and we have V* = - (lattice "divergence"). These geometric definitions may be extended 
to arbitrary lattices that need not be cubic or periodic [fill. In the minimal Coulomb 
gauge, the lattice divergence of A vanishes, (V*A) X = 0. 

At a minimum of the minimizing function we also have the condition d 2 Fu(s)/ds 2 \ s= o > 
0. From (49) we obtain 

d 2 Fu(s)/ds 2 = (2AT)- 1 X: ti{[u x+l -u x ]{d/dsWU) x , - (»U)l ti } }. (59) 

x,i 

With 

A x ,,(s) = A[(^UU), (60) 

the positivity condition reads 

d 2 Fu{s)/ds 2 \ s=0 = N- 1 ]T tr{[^ +i - u x ] dA x ^{s)/ds }| s=0 > 0. (61) 

x,i 

To make explicit the geometrical meaning this condition, and for future use, we define 
the lattice gauge covariant "derivative" D(A) by 

[D(A)w] Xji = -dA X!i (s)/ds\ s=0 . (62) 

It is the matrix that maps an infinitesimal gauge transformation u into the correspond- 
ing first-order change in the coordinates. (The minus sign is to be coherent with the 
continuum formula 9 A = A — D(A)u.) Like the lattice gradient V, it maps site variables 
into link variables, and we have D(l) — V. Its dual D(A)* maps link variables into site 
variables, and is defined by (D*(A)V, uj) = (V, D(A)uj), for arbitary link variables V and 
site variables uj. For the coordinates A defined above, D(A) is given explicitly by 

[D(A)u] x ,i = t b [D{A)u] b X!i = (1/2) [( U Xti oj x+i - u x U Xti ) - (h.c.)] traceless , (63) 
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where U = U(A). 

The positivity condition d 2 Fu(s)/ds 2 \ s= o > 0, which holds at absolute or relative 
minima of the minimizing function, thus reads 

£(V<,(Z^)<,>0forall^. (64) 

We define the lattice Faddeev-Popov matrix 

M{A) = D*(A)V, (65) 

which maps site variables into site variables. For transverse A, one may verify that this 
matrix is symmetric, 

D(A)*V = V*D(A) for V* A = 0. (66) 

We conclude that at a minimum of the minimizing function, this matrix is also non- 
negative 

M(A) > for A G A. (67) 

It has a trivial null space consisting of x-independent eigenvectors, Vu; = 0. This is a 
reflection of the fact that the minimal Coulomb gauge does not fix global (x-independent) 
gauge transformations. On the orthogonal subspace, M(A) is strictly positive for config- 
urations A that are interior points of the fundamental modular region A. (For additional 
properties of A, see fill .) 

10 Lattice Coulomb-gauge Hamiltonian 

In the preceeding section we obtained a parametrization of the gauge orbit space by 
means of transverse configurations restricted to the fundamental modular region. We 
shall now calculate the lattice Coulomb-gauge Hamiltonian H cou \ as the restriction of the 
Kogut-Susskind Hamiltonian H^s to gauge-invariant wave functions, Q( 9 U) = $(U). 

We represent states $(U) as functions of the coordinates, $ = $(A). The coordinates 
transform according to 9 A b x i = A b xi ( 9 U) under local gauge transformation g x . For an 
infinitesimal gauge transformation g x = exp(u x ) = 1 + u x , we have 

'<< = <i " [D(A)u]l ti , (68) 

by the definition (|6^) of the lattice gauge-covariant derivative D(A). With Q( 9 A) = 
$(A — D(A)w), we conclude that gauge-invariant wave functions, which satisfy $>( 9 A) = 
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<&(A) for all g, satisfy the first order differential equation, 



Y,(D(A)uy x ^/dA^ = 0, 



(69) 



for all lo. This is equivalent to the condition 



(D*(A)d/dA) b x & = 



(70) 



for all lattice sites x, which is the lattice version of Gauss's law constraint. This equation 
holds for any coordinates A Xti on the group manifold, provided only that the lattice 
gauge-covariant derivative D is defined by (|62|). 

We now proceed as in the continuum theory, and introduce the lattice Faddeev-Popov 
identity 



The ^-function 8& is defined by strict analogy to the continuum theory. The product 
Y\ x dg x extends over all sites x of the lattice, where dg x represents Haar measure. However 
the primed product, Ilx'^A[(V*A) ;r ] extends over all but one site xq of the lattice because 
the lattice divergences satisfy the identity J2x(^*A) x = 0, and are thus not all linearly 
independent. The primed matrix M'(A) is the lattice Faddeev-Popov matrix M(A) = 
D*(A)V, with the rows and columns labelled by xq deleted. One may show that the 
integrand is independent of xq and that, apart from an overall normalization constant, 
det M'(A) for V*A = is the determinant of M(A) on the space orthogonal to its trivial 
null space. In the minimal Coulomb gauge M'(A) is a strictly positive matrix, in the 
interior of the fundamental modular region A, and its determinant is positive. 

[The above Faddeev-Popov identity also holds for other coordinates, such as that 
provided by the exponential mapping U Xj i = exp A Xti , with gauge condition (V*A) X = 0. 
We now derive an explicit expression for the lattice gauge-covariant derivative D(9) in 
any coordinate system 6 xi , where D(9) is defined by the linear change in the coordinates 
induced by an infinitesimal local gauge transformation. We start with the Maurer-Cartan 
differential 



which relates infinitesimal changes in the coordinates d6 to infinitesimal changes dll in 
U. For an infinitesimal gauge transformation g x = 1 + uj x we have 




(71) 



U X: idU xi 



t a r a (9 Xii )de: 



(72) 



9 U X:i = U Xji + uJ x U X)i — U X)i uj x 
dU X: i =uj x U X: i — U x .iu x+ i 
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U x ,%dUl ti = U x ,iU x +-Ul ti - c^, (73) 

where we have used uj x = —uj x . The last expression represents the difference between the 
parallel-transport of u, from x + i to x, and uj x . By analogy with the continuum theory, 
it is natural to call this quantity the lattice gauge-covariant derivative in the site basis, 

[V{U)uU = U X ^ X ^U - ( 74 ) 

If we choose a basis in the Lie algebra, u x = t a u x , we have U Xi it?Uli = t a R ab {U), where 
R ab (U) is the adjoint representative of U. In the same Lie algebra basis we have the 
expansion [D(U)cu] X j = t a [V(U)u] x t . This gives the explicit expression for the lattice 
gauge-covariant derivative in the site basis 

PW<i = R a \U)u h x+l - u a x . (75) 

By analogy with the continuum formula for an infinitesimal gauge transformation 9 A = 
A — D(A)u, we define the lattice gauge-covariant derivative in the link basis Du as the 
first order change in the link coordinates induced by an infinitesimal gauge transformation 

dQl % ee [£>(*)<< • (76) 

The two quantities are related by use of the Maurer-Cartan differential 

t a r a (0 x ,)do:, = u x M,i = nmuu (77) 

which gives 

<(MW)<,= [^)<*- (78) 
The last formula may be inverted using Jf {0 Xt i)ip^{6 Xt i) = 6@, 

and we obtain for the lattice gauge-covariant derivative in the link basis the explicit 
formula, valid in any coordinate system 

[D(9Uli = J:(e xt )[R a \0K +i - "SI , (80) 

where we have written R ab (6) = R ab {U(6)).} 

We shall now derive the lattice Coulomb-gauge Hamiltonian H cou \ in any coordinate 
system with gauge condition Y^,x,i(A x ,i — A x _i). We insert the Faddeev-Popov iden- 
tity into the formula for the inner-product ($, Because the Haar measure d\J x ^ = 
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det ip(A Xti )dA Xt i is invariant under left and right group multiplication, it is invariant under 
local gauge transformation. Consequently we obtain, just as in the continuum theory, 

(<£>,*)= N [ Y[dA x ^(A x ,)]X 5[(V*A) X ] det M'(A)^(A)^(A), (81) 

JA x,i 

where N is the finite volume of the local gauge group, and we have written ip for det -0. 
The primed product YlJ extends over all sites but one, as explained in the preceding 
section, and correspondingly for M' . 

The gauge condition V*A = may be solved by a fourier decomposition of A Xji on 
a finite, periodic, cubic lattice, using longitudinal and transverse polarization vectors for 
finite k, and including k = or harmonic modes. This allow to write 

Ar,i,a = X^ )QI — (Vo" a ) X) i, (82) 
where (V*A]^) X = 0. Here A Tr includes both transverse and harmonic modes, A Tr = 

A tr + A h_ j ntegral 

over link variables which represents the inner product becomes an 
integral over fourier coefficients. Because of 5(V*A), only the transverse and harmonic 
modes survive. We express the result as 

(<5>,^)=N [ rfA Tr n^(^S)detM'(A Tr )<l>*(A Tr )^(A Tr ), (83) 

Jk x,i 

where dA Tv designates the interal over all transverse and harmonic modes. 

We wish to express the Hamiltonian Hks = T + V as an operator on the reduced 
space of gauge orbits parametrized by A Tr e A. The potential energy operator V acts in 
the larger space by multiplication by a gauge-invariant function V( 9 A) = V(A). On the 
reduced space it acts simply by multiplication by V = V(A Tt ). 

The Lie generator J a is symmetric with respect to Haar measure, so we may write the 
kinetic energy operator T in the Kogut-Susskind representation as the quadratic form 

($,Ttt) = / n^G?o72) E(£m,«*)* ( 84 ) 

Only first derivatives of the wave-function appear here, which will allow us to apply 
Gauss's law as in the continuum theory. Here we have written 

£x,i,a = Jx,i,a (^^) 

in order to emphasize the analogy with the continuum electric field. The operator £ x ^ a 
represents the electric field in the the site basis. We also introduce the electric field in 
the link basis. 

E x ^ a = id/dA",. (86) 
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The two are related by 

£-x,i,a J a (A x> i) E Xj i :a . (87) 

In this notation, Gauss's law constraint reads 

{D{A)*E) X <& = 0. (88) 

The operator £ x ^ a = J x ,i, a transforms homogeneously under gauge transformation, so 
the integrand of the quadratic form T is gauge invariant. We insert the lattice Faddeev- 
Popov identity into this integral, and obtain 

(<t>,T*) = iV / HdA x ^(A x>i )detM'(A) 

J A 

x,i 

:ir*[(VM)J(flg/2) Efe,^)* ( 89 ) 



x 



We wish to use the Gauss's law constraint to express £ as a derivative operator that 
acts only on transverse variables. For this purpose we make a Fourier decomposition of 
the electric field E x ^ a on the periodic cubic lattice. For the A" -field we write the fourier 
decomposition, as 

A x = Y,<I>x,kAk, (90) 

K 

where X represents position and vector indices, and K represents wave vector and po- 
larization indices. Here the <f>x,K form a real orthonormal basis, and the inversion reads 

A k = Y / <Px,kA x . (91) 

x 

We have 

d ^dA K d ^ d 

dT x -\dr X M- K -\^ K ^ (92) 

so the fourier decompositions of A x and of J^- are the same. For the electric field E x 
we make the decomposition onto the same basis 

E x = Y j <Px,kE k . (93) 

K 

Consequently, from 

£ * = W < 94) 

we obtain 

E K = zL-. (95) 
dA K 



17 



We choose a polarization basis which is adapted to the decomposition into transverse and 
longitudinal and harmonic parts. Corresponding to the decomposition A = A Tr — Vex, 
we have 

E = E Tt - Vfi, (96) 

where V*A Tr = V*E Tr = 0, and we have introduced the lattice color-Coulomb potential 
operator fi. It follows that if A Tv has the fourier decomposition 

4* = E <Px,kA k , (97) 

K<=T 

then E Tv has the fourier decomposition 

El T = tJ2 (98) 

KeT OA K 

where T is the set of transverse and harmonic modes K. This expression defines E Tr as 
a differential operator that acts within the space of functions of A Tr . 

We shall now use the Gauss's law constraint (D*E) X Q = to determine the color- 
Coulomb potential Q x & on the reduced subspace, parametrized by A Tr e A. We write 
Gauss's law in the form 

D*(E Tr — VQ)<E> = 0, (99) 

which gives 

zrvfi$ = p$, (100) 

where 

(p) x = (D*E Tt ) x = [(D - VyE T % . (101) 

This operator represents the color-charge density carried by the dynamical degrees of 
freedom of the gluons. We shall only need to solve Gauss's law for fi$ in the reduced 
space where the Faddeev-Popov matrix is symmetric in the minimal Coulomb gauge, 
M(A Tr ) = D*(A Tr )V = V*D(A Tr ). Thus we have also 

(V*D Q.) x $ = Px $. (102) 

Because Y^x(^*V)x = for any link variable V, the last equation is consistent only if the 
total charge Q = J2 X Px of the state $ vanishes. Thus for gauge-invariant states in the 
minimal Coulomb gauge we have the constraint on the reduced space 



Q$ = 0. (103) 
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In the preceding section we have seen that in the lattice Coulomb gauge D*(A Tr )'V is 
also a strictly positive matrix for A Tr in the interior of the fundamental modular region 
A. Consequently the equation for Q always has a solution, 

a $ = (zrvrV $ . (104) 

With p = D*(A Tr )E Tt , this expresses the color-Coulomb potential operator Q on the 
reduced space as a derivative with respect to the components of A Tr . We conclude that 
in the reduced space the color-electric field E = E Tr — Vf2 acts according to 

where E Tr is defined by its fourier expansion ( |98|) . 

With this result we may express the kinetic energy operator T as a quadratic form 
on the reduced space, 

($, TV) = N f dA Tr J] ^{A^) det M'(A Tr ) 

Ja x,i 

x(g 2 /2)J2(£^Y (106) 

where $ = $(v4 Tr ) and similarly for \I>. Here £ $ is defined by £ x ^ a $ = J^A^AE,,^^, 
and E§ is defined by the preceding equation. This completes the specification of the 
lattice Coulomb-gauge Hamiltonian, H cou \ = T + V. 



11 Conclusion 

We have derived by elementary methods the lattice Coulomb-gauge Hamilonian H con i 
from the Kogut-Susskind Hamiltonian i^KS- The essential step was to write ifxs as a 
quadratic form, involving only first derivatives, so Gauss's law D*E<& = could be use 
directly to solve for the color-Coulomb potential. 

We have not addressed here invariance with respect to finite gauge transformations, 
which involves identification of points of the boundary of the fundamental region A, so 



that Gribov copies are avoided. These questions are addressed in |[L1|| , and have an 
important influence on the spectrum fl], ||, ||, [§] . 

However we believe that Gauss's law is essential to the phenomenon of confinement 
in QCD ||. We have shown how, in Hamiltonian theory with lattice regularization, 
the imposition gauge invariance in the form of Gauss's law leads to an instantaneous 
color-Coulomb potential in the lattice Coulomb-gauge Hamiltonian H cou \. 
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